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Abstract 

The present paper is devoted to studying the Hill-type formula and Krein-type trace formula 
for ODE, which is a continuous work of our previous work for Hamiltonian systems [[4] . Hill-type 
formula and Krein-type trace formula are given by Hill at 1877 and Krein in 1950’s separately. 
Recently, we find that there is a closed relationship between them [¥]. In this paper, we will 
obtain the Hill-type formula for the ^-periodic orbits of the first order ODEs. Such a kind of 
orbits is considered naturally to study the symmetric periodic and quasi-periodic solutions. By 
some similar idea in [4], based on the Hill-type formula, we will build up the Krein-type trace 
formula for the first order ODEs, which can be seen as a non-self-adjoint version of the case of 
Hamiltonian system. 
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1 Introduction 

In the present paper, we will study the Hill-type formula and Krein-type trace formula for S- 
periodic solutions for the first order ODE. Hill-type formula was introduced by Hill [2] when he 
considered the motion of lunar perigee at 1877, and the Krein-type trace formula was built up 
by Krein is m in 1950’s when he studied the stability of Hamiltonian systems. Although they 
appeared separately, there is a closed relationship between them. In fact, the Krein-type trace 
formula could be derived by the Hill-type formula for 5'-periodic orbits in Hamiltonian systems. 
Moreover, motivated by the Krein’s original work, the Krein-type trace formula was used to study 
the stability problem in n-body problem, details could be found in [am- In the case of Hamiltonian 
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system, the corresponding differential operators are self-adjoint, and the formulas for the first order 
ODE can be seen as a non-self-adjoint version. 

Let M m ( C) be the set of m x m matrices on C m , and denote by 

<B{rn) = C([0,T\,M m (C)), (1.1) 

the set of continuous path of m x m matrices on [0,T], We consider the n-dimensional first order 
ODE with the S-periodic boundary problem 

x(t ) = D(t)x{t), (1-2) 

®(0) = Sx(T), (1.3) 


where S is an orthogonal nxn matrix and D E 53(n). It is natural to study the S'-periodic solutions 
of the first order ODE when we study the symmetric periodic solutions and quasi-periodic solutions. 
In what follows, we always denote by 7 £>{t) the fundamental solution of the first order ODE (11.21) . 
that is, 7 nit) = D(t)^ D {t) with 7 d(0 ) = I n - 

Consider ^ as the unbounded closed operator densely defined on L 2 (0,T;M n ) with the domain 
D S = { z(t ) € W 1 ’ 2 ([0,71,C") | z(0) = Sz(T)}. 

D is a bounded operator acting on L 2 (0,T;M n ) defined by ( Dz)(t) = D(t)z(t). In this paper, we 
will prove the following Hill-type formula. 


Theorem 1.1. Let D € 53 (n) and S be an orthogonal matrix. Then, for any i/£C 


det 



D + vl r 


dt 


+ Pc 


-In 


(-l) n |C(S)|e"^e-^ io T T < D ) dt det(S 7 d (T) - e uT I n ), (1.4) 


where Po is the orthogonal projection onto ker(S — I n ) and C(S ) is a constant depending only on 
S. 


Here, if we let W = ker(S — In) 1 - and ko = dim ker(S — I n ), then C(S ) = T k ° 


_ 1 . 

det[(5— I n )\w] ' 


Remark 1.2. 1) In Theorem |T7Tj since {DPq — nl n )y^ + PqJ is not a trace class operator, 

but a Hilbert-Schmidt operator. Hence, the infinite determinant 


-1 


det 


W d 


\dt 


— D + v I r 


+ Po 


-1 


= det 


id-(D + P 0 -nI n )(j t +P 0 ' 


is not the classical Fredholm determinant. In fact, we can define the determinant in the 
following way. Let Pn be the orthogonal projections onto 


V N = 


© 


I 'ecr(^),M<A r 

And the conditional Fredholm determinant 




det 


r/ d 


Adt 


— D + ul r 


A Pi 


-ii 


= lim det 
iV—>00 


id — Pn(D + -Po ~ vln) \-^ A Po 


-1 


N 
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2) If S = I n , then the boundary condition problem ll.Hil.3 j) is the canonical periodic boundary 
condition problem. In this case C(S ) = T~ n . When the period T = 1, the Hill-type formula 
was obtained by Denk W- Based on the Hill-tpye formula, Denk developed some efficient 
numerical method for the ODEs. 

Similar to (4j, the Hill-type formula (11.41) is the starting point of the Krein-type trace formula. 
To get the trace formula, for Do(t), D(t) € 03 (n), we consider the eigenvalue problem 

z(t) = (D 0 (t) + aD(t))z(t), (1.5) 

z{ 0) = Sz(T), (1.6) 

that is, to find the a € C such that the system (|1,5111.61) has a nontrivial solution. As above, let 7 a 
be the fundamental solution of (11.51) . To state the trace formula, we need some notations. Write 
M = 5yo(T) and D(t) = 7 c j' 1 (t).D(f) 7 o(t). For k £ N, let 

M k = f D(ti) [ D(t 2 )••• [ D(tk)dt k ■ ■ ■ dt 2 dt-i, 

Jo Jo Jo 

and 


G k = M k M (M - e vT I 2n ) 1 . 

Theorem 1.3. Let fgC such that ^ — Dq — v is invertible, F = D (^ — Dq + vl n ) 1 , then 

Tr{F) = \ [ T Tr(D(t))dt - Tr(Gi), (1.7) 

1 Jo 


and for any positive integer m > 2, 


Tr(F m ) - 


(- 1 )' 


k =1 


k 


E T Va n ---a lk ) 

ilH-1 -jk=m 


( 1 . 8 ) 


Remark 1.4. (1). For m = 1, F is not a trace class operator but a Hilbert-Schmidt operator. 

And hence Tr{F) is not the usual trace but a kind of conditional tracef^. That is 


Tr(F) = lim TrPrsrFPr^. 

N—yoo 


(1.9) 


For m > 2, F m are trace class operators. Obviously, for v = 0, Aj is the eigenvalues of 
U.5\ )- U.bV . if and only if j- is an eigenvalue of F. 


E 


1 



E 


k =1 


k 


E Tr(G h -G jk ) 

j lH-1 -jk=m 


where the sum takes for eigenvalues counting algebraic multiplicity. 


(2). The trace formula for ODE is different from that for Hamiltonian system, and the reason 
is that the differential operator here is not self-adjoint any more. 
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Remark 1.5. The idea and the techniques are similar to those in however, there are at least 
two reasons to write this paper. For the first, it is not like that Hamiltonian system comes from 
mechanic system mostly, the ODE systems come from many areas, then the Hill-type trace formula 
and Krein-type trace formula for ODE will be more convenient to be used. For the second, since 
the differential operators for ODEs are not self-adjoint, we should do some spectral analysis for 
ODE carefully. Although the formulas for ODE are similar to that for Hamiltonian system in B 
however, it is difficult to deduce them. 

The trace formula for Hamiltonian systems is a useful tool in study the stability of Hamiltonian 
systems, by using the trace formula and Maslov-type index theory [8], some applications to the 
n-body problem is given in W In Section 5, as an application, we will give a generalization of 
Krein’s work on second order systems [6j. 

This paper is organized as follows, in section 2, we review the basic properties of conditional 
Fredholm determinant and conditional trace. In section 3, we derive the Hill-type formula for the 
S'-periodic orbits in ODE. In section 4, we get the trace formula from the Hill-type formula. Finally, 
as an example, we will reformulate Krein’s trace formula from our viewpoint. 

2 Preliminaries 

In this section, we will mainly recall some fundamental properties of conditional Fredholm deter¬ 
minant, which was developed in [5]- In the classical settings, the Fredholm determinant det(id + -F) 
is defined for a trace class operator F, details could be found in [TO] , However, when we study the 
ODEs, the operators we encountered are of the form (id + F) with that F is not a trace class op¬ 
erator, but a Hilbert-Schmidt operator. Thus, the Fredholm determinant det (id + F) could not be 
defined, and the conditional Fredholm determinant will be used instead. To define the conditional 
Fredholm determinant , the trace finite condition plays an important role. 

Let {Pk} be a sequence of finite rank projections, such that the following conditions are satisfied, 

(1) for k < m, Range(Pk) C Range(P m ), 

(2) Pk converges to id in the strong operator topology. 

A Hilbert-Schmidt operator F is called to have the trace finite condition with respect to Pk, if the 

limit lim Tr(PkFPk) exists, and the limit is finite. Clearly, a trace class operator has the trace 
k—yoo 

finite condition. Obviously, all the Hilbert-Schmidt operators with trace finite conditions consists 
a linear space, that is, if F\ and Fi are Hilbert-Schmidt operators with trace finite condition, then 
a\F\ + (X 2 F 2 has the trace finite condition. 

By [10], if F is a Hilbert-Schmidt operator, then the regularized Fredholm determinant is defined 
by 

det 2 (id + F) = det id + F)e~ p ^j . (2-1) 

As been pointed in [5], if F is a Hilbert-Schmidt operator with trace finite condition, then the 
conditional Fredholm determinant can be defined by 

det (id + F) = lim det (id + PkFPk ) 

k —^00 

= det 2 (id + F) lim e Tr ^ PkFPk \ (2.2) 
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where, P k FP k are finite rank operators, and hence det (id + P k FP k ) is well defined. As we proved 
in [5] , many fundamental properties of conditional Fredholm determinant are similar to that of the 
usual Fredholm determinant. 

Proposition 2.1. 1) If F\ and F 2 are Hilbert-Schmidt operators with trace finite condition, then 

det ((id + F\)(id + F 2 )) = det (id + F\) det (id + F 2 ). 


2) Let E = E\ © E 2 , and F{ be Hilbert-Schmidt operators on Ei with trace finite condition with 
(i) 

respect to P k , i = 1,2. Let F = F± © F 2 , then F has the trace finite condition with respect 
to P^ © Pjf\ and 

det (id + F) = det^c^ + Fi) det (idE 2 + F 2 ), 
where idp, are identities on Ei, for i = 1,2. 

Similar to that we have given in [5], it is not hard to show that det(?'d + aF) is analytic on a, 
for a Hilbert-Schmidt operator F with trace finite condition. For reader’s convenience, we will give 
the proof of it. 

Lemma 2.2. Let F be a Hilbert-Schmidt operator with trace finite condition with respect to {Pk}, 
then det (id + aF) is an entire function. 

Proof. Write 


f k (a) = det (id + aP k FP k ), 


(2.3) 


then, by the definition, fk converges to f(a ) = det(*d + aPkFPk ) point-wisely. Moreover, since 
P k are finite rank projections, P k FP k are finite rank operators, and hence they are trace class. 
It follows that fk(ct) are entire functions. By Montel’s Theorem, it suffices to show that, {fk} is 
locally bounded, that is, for any compact set K C C, there is a constant C > 0 depending only on 
K such that 

sup{|/ fc (a)| | a € K) < C. 


By 13 Lemma 2.3], 


f k (a) = det 2 (id + aP k FP k )e- aTrP * FP K 


Firstly, by [TO] Theorem 9.2], 

1 n „„ mo C||F||? sup H 

sup |det 2 (id + aP k FP k )\ < sup e c W aPkFPk H* < e “G* 

aeK a&K 


where II - M 2 is the Hilbert-Schmidt norm. Secondly, since lim TrP^FP^ exists, we have that, 

N —>00 

there is a constant C\ > 0 such that for any k € N, \TrP k FP k \ < C\. It follows that 


sup |e aTrPkFPk 

OidK 


Ci sup |a| 

< e aeK 


Therefore, 

sup {|/fc(a)| | a G K) < C 


for some constant C which depends only on K, that is, {f k } is locally bounded. The proof is 
complete. □ 
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Remark 2.3. In the proof of Lemma \2J\. we show that {det (id+aP k FP k )} is a normal family, and 
hence there is a subsequence of {det(id + aP k FP k )}, say {det (id + aP kj F P kj )}, which is convergent 
uniformly to det (id + aF) on any compact subset of C. Denote by 

gkj (a) = det (id + aP kj FP kj ). 

Following m , we call a function / : X —» Y between Banach spaces, finitely analytic if and 
only if, for all A\, ...,A n € X, f[z\A\ + ... + z n A n ) is an entire function of z \,..., z n from C n to 
Y. This concept is very useful in studying the conditional Fredholm determinant. An important 
property is the following: 

Theorem 2.4. 1 10. pp-45, Theorem5.1J. If a finitely analytic function f satisfied f(x) < G(|| x ||) 
for some monotone function G on [0, oo), then f is Frechet differentiable for all x € X, and Df is 
finitely analytic function from X to £(X,Y) (Banach space of the linear operators from X to Y), 
and 

(Df)(x)<G(\\x\\ +1). 

Next, we will consider the Taylor expansion of det (id + aF). Write 


g(a) = det (id + aF). 

By m Theorem 5.4], we have the following lemma. 

Lemma 2.5. Let g k (a) = det (id + aF k ). Then the Taylor expansion near 0 for g k {a) is 

OO 

g k (a) = ^2 a m a k , m /m\, 


(2.4) 


where 


det 


m =0 

( TrF k m — 1 0 

Tr(F^) TrF k m — 2 

Tr(F™~ l ) Tr(F™~ 2 ) ••• 

V Tr(Fff) TffFff- 1 ) ... 


0 \ 
0 

TrF k 1 
Tr(F 2 ) TrF k 


and F k = P k FP k . 


The following reasoning is similar to that in [U Section 2], and we write here again for reader’s 
convenience. 

Now, let g(a) = X) TJ 0 ™ be the Taylor expansion of g(a). Since g kj (a) converges to g{a) on 

m 

any compact subset of C, the coefficients a kjt m —>• a m . Notice that F is a Hilbert-Schmidt operator 
with trace finite condition, then TrF k —» TrF. Therefore 

/ TrF m — 1 • • • 0 \ 

Tr(F 2 ) TrF ■■■ 0 


a m — det 


\Tr(F m ) Tr{F m ~ l ) ■■■ TrF 
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Note that for a small, by m p.47,(5.12)], we have that 

det(zd + aFk) = exp ^ ---- a n, Tr{F™') \ . (2.5) 

m= 1 

By taking limit, we have the following theorem. 

Theorem 2.6. Let g(a) = det (id + aF). Then for a small, 

g{a) = exp ( ]T - a m Tr{F m )) . (2.6) 

m= 1 

At last of this section, we will brief review the Hill-type formula for Hamiltonian systems. 
Assume 5 be a symplectic orthogonal matrix, for 5-periodic orbits in Hamiltonian system, the 
Hill-type formula was proved in [5], which is listed as follows. 

Theorem 2.7. There is a constant C(S) > 0, which depends only on S, such that for any v € C 

^ (( ~ J Jt ~ B ~ UJ ) ( _ J Jt + P0 ) _l ) = C ^) e ~^ J ° Tr{JB{t))dtX ~ ndet (^ T ) - A/ 2n), ( 2 ‘ 7 ) 

where A = e uT and C(S) > 0 is a constant depending only on S. 

Precisely, 


n 

C{S) = 2~ n (2/T 2 ) k ° [| 

j=k o+l 


1 

1 — cosh(-v/—1 VjT) ’ 


where fco = dim ker(i? — yj—IQ — I n ). In fact, if we denote W = ker(5 — hn)' 1 , then 


C(S) = T~ 2k ° 


1 

det(S — l2n)\w 


Obviously, C(S) >0. If S = l 2 n , he. for the periodic boundary conditions, C(S) = T 2n , and if 
ker(S — l 2 n) = 0, then 


C(S) = det (5 — hn) 1 - 

For a n x n orthogonal matrix S, let W = ker(S' — /„)-*-, and we also denote 

C(S) = T~ k °- 


1 


det(5 — I n )\w ’ 

where k q = dim ker(S’ — I n ). It can be shown that C(S) = (—l) 0- 10(5)1, where a is the orientation 

'so' 


of 5. Then in the special case 5 = 


0 n 5 


, where 5 is a n x n orthogonal matrix, 


C(5) = 0(5) 2 . 


( 2 . 8 ) 


3 Hill-type formula for ^-periodic solutions of ODEs 

In this section, we will mainly derive the Hill-type formula for the first order ODE. 
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3.1 Hill-type formula for the first order ODEs 

To derive the Hill-type formula, we will consider the conditional Fredholm determinant of ODEs. 
Firstly, we consider the first order ODE with the S-periodic boundary condition, 


z(t) = D(t)z(t ) (3.1) 

*(0) = Sz(T), (3.2) 


where S is an orthogonal matrix on M n and D € 95 (n). To continue, we will consider the spectrum 
of with the S'-boundary condition. Denote the domain of ^ by 

D S = {z(t) € W 1,2 ([0,T],C n ) | *(0) = Sz(T)}. 

Obviously yf—1 ^ is a self-adjoint operator on E = L 2 ([0, T], C 2n ) with domain Ds- 

For simplicity, write A = \J — 1 (jj . By some simple calculation, we have that v € cr(A) if and 
only if ker (Se'^~^ vIriT — I n ) is nontrivial, and the corresponding eigenvector is where 

£ € ker {Se^~^ vInT — I n ). Since S is an orthogonal matrix, there is a unitary matrix U such that 


/ e -V=I*i 


u*su = 


O— V~ 1®2 


0 -V=ie n 


where 0 = 6\ = • • • = 6 ko < 9 ko + 1 < ■ ■ • < 0 n < 27r. It is easy to check that v € ct(H) if and only 
if v = (0j + 2kiT)/T, for some 1 < j < n; furthermore, (U*SUe^~^ u ^ InT — I n )f,j = 0 if and only if 
= (0,0, • • • , 0,1,0, • • • ,0). We have the following lemma. 


i -1 


Lemma 3.1. The spectrum of ^ (with the S-boundary condition) is periodic with the period of 
2n/T. Precisely, let Vj = 0j/T, then 


7 

cr(—) = IJ { - \/-l {vj + 2/c7t/T) | k € Z}. 
i=i 


/n /act i/j = Qj/T. Moreover, let — + 2kTr/T) € cr(^), f/ien f/ie corresponding eigenvector 

is e (^+2W where = (0,0, - - - , 0,1,0, - - - , 0). 


j-i 


It is worth being pointed out that 


0 = Vi = ■ ■ ■ = u ko < u ko+ i < • • • < v n < 2n/T, 


(3.3) 


where ko = dim ker (S’ — I n ). 


Remark 3.2. Since S is an orthogonal matrix, e ie is an eigenvalue of S if and only if e~ l9 is an 
eigenvalue of S. By the argument before Lemma [37l\ we know that, under the S-boundary condition, 
the spectrum 


a 




(3.4) 


which will be used later. 




-1 


By Lemma 13.11 + PqJ is a Hilbert-Schmidt operator, where Po is the orthogonal projec¬ 

tions onto ker (^). Applying a similar reasoning to [5J Remark 2.9] shows that (D + Po)(^ + Po) _1 
has the trace finite condition with respect to {Pat}, where P/v are orthogonal projections onto 


Vn= © ker ( F -5 


(3.5) 


UGa (di)’\ U \< N 

In fact, we have the following lemma. For reader’s convenience, we will give the proof of it. 
Lemma 3.3. Under the above assumption, we have that 


d 1 ^ f T ^ 1 + cos T Vi f T 


Pi Tt ( PnD {i; +P o ) lp ") = rP[ Uidt+v^i Y 


dt 


T ~r'« 


j=k 0 +1 


2 sin T v 


3 Jo 


Djj{t)dt, 


where D = U*DU. 

Proof. For 1 < j < n, let 


J K. = ® tar (|+' /z + 

kez 


(3.6) 


then E = ®Mj. Let -Le(‘ , h 2i ” r / T )^" t !7^ be the eigenvector for —^ with respect to the 
eigenvalue —a/—1 (vj + Then 


lim Tr(P N D(^ + P Q ) ' P N 

N=>oo \ \dt J 


V lim V (d(^- + P { 

L—d N=>-oo \ \dt 

3=1 |fc|<iV X 


^ ^ 1 \ p {yj+2kir/T)y/^Il n tjj£' _}_ p (uj+2kir/T)P^lInt 


01 7F 


•3 > 






For 1 < j < ko, 




For fco + 1 < J < n, 

lim V (d(P- + P { 
N—>oo \ V /7/ 


|/c|<AT 


^ 1 ^ (v,+2fc7r /T)y/^\I n tjjt 1 (vj+2kir/T)y/^\I n t 


dt ' * Vr e 


‘^•7T e 


C/{, 


2 lim ------— [ 

N ^°° UPij V=Hyj + 2fcvr/T) J 0 


k<N v ■? 

.—-1 + cos TVj f T 
= v —1 1 


+ 2kir/T) Jo 
Djj(t)dt. 


Djj(t)dt 


sin T i/j Jq 

The above calculations imply the desired result 


□ 
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Now, we will consider the conditional Fredholm determinant det 
respect to Pn . As we have done before, write 


&-D)(i + P 0 ) 




with 


— D 


+ Po 


-l 


= id - (D + P 0 


+ Po 


-l 


Hence, the conditional Fredholm determinant det I id — (D + -fo)(4 + Pq) 1 I is well-defined 


Please note that for any D\ E 25(n) such that 4 + D\ is invertible, the operator D(4r + D\) 1 isa 


dt 1 10 ^\dt 

Hilbert-Schmidt operator with the trace hnite condition with respect to Pjy. Therefore the infinite 

-ll 

is well defined. 


determinant det 


i-D 


d 


In this remaining part of this subsection, we will deduce the Hill-type formula for first order 
ODE with S-periodic boundary condition, where S is an orthogonal matrix. Consider the following 
equation 


u(t) = D(t)u(t), 
u(0) = Su(T). 

To obtain the Hill-type formula of the above system, we will consider the following Hamiltonian 
system, 


where B{t) = V 


iD 


V~\ S = 


z(t ) = JB(t)z(t), 
z{ 0) = Sz(T), 

' S 


S 


. Changing the basis by V, we have 


(3.7) 

(3.8) 


d 


V * [ - J dt ]V = 


— l dt 


—i 


■ d ),V*B(t)V = 


iD 


dt 


Or, 


and V*JV = 


— Hr, 


iPi 


Let 7 (t) be the fundamental solution of (13.71) . Under the new basis, it is obvious that 

7£>(*) \ 


7 (f) = 


P 


where 7 u(t) satisfied 725 (f) = with 725 ( 0 ) = p. 

By the Hill-type formula (|2.7I) for Hamiltonian system, we have 


det 


l dt 


—1 


• d 


dt 


= C(S)e nvT exp 


iD 


Tr 


—iu I„ 




iuP 


-iTt+ P o 


D 


0 , 


| dt 


det (5 7 (T) - e vT l 2n ) 


By Proposition 12.1( 2). we can rewrite the above equation as 


det 


d 


i— - iD + iuP ) ( i— + Po 


dt 

= C(S)e~ nvT exp 


d 


det 


d \ d 

-p-™ 1 " \- i dt + p “ 


-1 


Tr{D)dt 


det (S 7 25 (T) - e uT p) det (S - e vT P) . (3.9) 
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(4 - iD + (4 + A) 


-i 


-i 


To calculate the determinant det 
terminant det (—4 — (“4 + 4)) 

R n_fc °, obviously 5 = | ^ j . We only need to compute det 


, it suffices to calculate the de- 


. Let Si = Slw, which is an orthogonal matrix on 


(-4 - ivI ko) (-4 + £>) 


-1 


with T-periodic boundary condition and det 


■ dW 1 


( *dt wln-ko) ( 4) 

A 

dt 

d \ f 2kir 


with Si boundary condi¬ 


tion respectively. Firstly, notice that the spectrum of i4\ with T-periodic boundary condition 


a l dt 


X T 


-,k € Z > . 


Moreover, by direct computation 


n< i + 


IV 


2fc7r/T ' 


= iv 




fceN 


(2fcvr/T) 2 ' 


2 * 

T 

i 


sinh(Ty/2) 


= - 1). 


We have 

det 


(-iL - ;„ 4o ) (-;L + P 0 ) ‘ (3,10) 


Secondly, writing Si = 


Si 


Si 


|, by the Hill-type formula for Hamiltonian system with Si 


boundary condition, we have 
det 


4 + ^ 


-4-*" 


iA 

L dt 


-1 


Therefore 

det 


l dt 


= C(S 1 )e-( n - k0 > T det (Si - e vT I n - k ,)‘ 


^ dt 


) ( _ *^) =vC(Si) 2 e ( 2 0) det (Si — e uT I n _ ko ) , (3.11) 

where <7 can be ±1, to be decided. Please note that, when we take v = 0, the left hand side of 
(13.lip equals to 1 and the right hand side equals to erdet (Si — I n -ko) C'(S) 2 . Since C(S) > 0 and 
a is a square root of 1, we have 

cr = sign( det (Si - I n -k 0 ))- 

Since C{S\) = [det ((Si — 2 , it follows that aC{S\)^ = det ((Si — I n - ko )~ l ), and hence 


det 


i d j i'L'In—ko 


-l- 


, d 
dt 


( n—kr\)vT 

= C(Si)e-det (S - e" T / n _ fc0 ) , 


(3.12) 


where C(Si) = det ((Si — I n - ko ) 4 From (13.101) and (13.121) . we have the following lemma. 
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det 


Lemma 3.4. Let S be an orthogonal matrix. Then 

4 - iw ”) (~P + p °) 

where C(S) = 

Substituting (13.131) in (13.91) . we have the following proposition. 

Proposition 3.5. Under the above condition, we have 

-T 


= (-1 ) k H ko C{S)e~ s ^- det (5 - e vT I n ) , 


(3.13) 


det 


d 


i— -iD + wl n i— + Pt 


dt 

Notice that 


d 


dt 


= {-l) k °i ko C{S)e-^e-^^ Tr ^ D)dt &et{S lD {T) -e uT I n ). 


det 


d 


i— — iD P wl n i— + P( 


dt 


. d 


= det 


= det 


= det 


dt 


— Dp vl r 


— D P vl n ) ( ~ + Po 


dt 

d_ 

dt 

d 


-i 


- iP 0 


-l 


dt 


— Dp vl r 


+ Pc 


det 


Po 


dt 


~77 + Pq ) ( ~77 ~ iP0 


dt 


-l 


where the second equality holds true because of Proposition ^. If 11. Therefore, we have the following 
theorem which is just Theorem ll.il 


Theorem 3.6. Under the above assumption, 


det 


d 


~dt. ~ D P v ^ n 1 > ~ + P( 


d 


dt 

ivT 1 rT , 


-1 


Please note that 


thus we have 


= (-l) k °C(S)e-^e- 2 fo Tr W dt det(Sp D {T) -e vT I n ). 
(-1 ) k °C(S) = (-l) Tl |C'(S’)|, 


(3.14) 


det 


d \ / d 

--DPuI n j (-PP 0 

nvT 1 r T 


-1 


= (-l) n |C'(5)|e- I ^-e-2/o Tr (D)dt d et(5 7D (r) - e uT I n ). 


det ( ( -J&-B l-Jf t + P 0 


Remark 3.7. From f5j, ' ' jd 
function G. By the above analysis, we have 

-l 


-l 


< G(|| B ||) for some monotone 


det 


— Dp vl r 


+ Pc 


= det ( 

-4-b)( 

V 

v dt J \ 


dt 


-i 


• d(v), 
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where d{v) is a function depending only on v and B(t) = V 


iD 


V 1 . Thus we have that 


det 


(a-0(4 + ^) 


-1 


is a finitely analytic function from 03 (n) to C, and satisfied 


det 


d \ (d A 
dt ~ D ) [dt + 0 


-l 


< Gi(|| ||), 


for some monotone function G\. Moreover, from Theorem \2.f\ suppose Ll C C m be an open set, 
and D(Z) is analytic map from Ll —»• 03(n) ; f/ien det ~ D(Z)^ + Po^j ^ is an analytic 

function on 

By the multiplicative property of conditional Fredholm determinant, Proposition 2.1(1) , we 
have the following corollary. 


Corollary 3.8. Let D € 05(n) such that ^ — D is invertible, then 


det 


dt 


— D\ + vl r 


T-o" 1 

dt 


= e - nuT / 2 e -i So TpD-D^dt det (s 1d (T) - e vT I n ) det (S 7Dl - 4)" 1 • 


(3.15) 


4 Trace formula for 1st order ODE 

In this section, we will derive the trace formula from the Hill-type formula. Firstly, we will study 
the Taylor expansion for linear parameterized Monodromy matrices, which is similar to the case of 
Hamiltonian systems [3]. 

4.1 Taylor expansion for linearly parameterized Monodromy matrices 

Let Dq,D € 05(n). For a € C, set D a = Dq + aD , for a € C, let 7a be the corresponding 
fundamental solutions, that is 


fait) = D a (t)j a (t). 

Fixed no € C, direct computation shows that 

a(t)) = ^{t)(D a (t)~ D ao (t))^ a (t) 

= ( a - oc 0 )^(t)D(t)'y ao (tyy^(t)'y a (t). 

Without loss of generality, assume ao = 0. In what follows, write 

7 oft) = 7o ' 1 ^ba(4), 


and 


D(t) = 7o 1 (t)H(t) 7o (t), 
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thus 


^7 a(t) = aD{t)%{t). (4.1) 

To simplify the notation, we use to denote the k -th derivative on a. Taking derivative on a 
for both sides of (14.11) . we get 

= D(t)%(t) + aD a (t)j^(t). (4.2) 

By taking a = 0, 7o(t) = I n , we have 

7?V) = [ D(s)ds. 

Jo 

Now, taking derivative on a for both sides of (14.21) . we get 

^7i 2) W = 2 D(t)^\t) + aD(t)^\t). 

Take a = 0, and we get 



£>(s)^ ) (s)da. 


By induction, 


^7o k) (t) = kD(t)*/k k l \t), 

7?V) = k [ D(s) A /Q k l \s)ds. 
Jo 


For t = T, by Taylor’s formula, 


%(T) = I 2n + a^\T) + • • • + a k ^ k) (T)/k\ + • • • , 


where 


7 o\T)= D(t)dt 


and 

i k) (T)/k\ = f D{t)^- l \t)/{k - 1 )\dt,k€ N. 
Jo 


By induction, we have 


# (T)/k\ = 


f T D( tl ) f 

Jo Jo 





D(tk)dtk ■ ■ • dt 2 dti, k € N. 


Obviously 7 Q (T) is an entire function on the variable a. We summarize the above reasoning as the 
following proposition. 
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Proposition 4.1. Let D a = Do + aD, j a (T) be the corresponding fundamental solutions. Write 
7 a = % lr )a- Then, the Taylor expansion for j a (T) at 0 is 

7 a(T) = In + a^\T) + • • • + a k $\T)/k\ + • • • , (4.3) 


where 


7 o k \ T )/k'-= f D(t i) f D(t 2 )••• [ D(t k )dt k - ■ ■ dt 2 dti,k 
Jo Jo Jo 

In what follows, to simplify the notation, set 


€ N. 


(4.4) 


M(a)=%(T), M 0 = In and M j = $\T)/j\, j G N, 


then 

OO 

M(a) = J2 aj Mj- 

3=0 

Set M = S' 7 o(T), then 57 Q (T) = MAI (a). For A E C, which is not an eigenvalue of M, by 
some easy computations, we have that 

det(5 7a (r) — AI 2n ) = det(MM(a) - A/ 2n ) 

= det(M - A I n + aMMi + • • • + a k MM k + • • •) 

= det(M - A I n ) det(/ n + • • • + a k {M - A I n )~ l MAI k + ■■■). 

Let 


G k = (.M - A I n )~ l MM k , 


(4.5) 


and 

/ (a) = det(/ n + • • • + a k G k + ■■■), 

which is an analytic function on C. Next, we will compute the Taylor expansion for /(a). Let 

OO 

G(a) = ak ~ 1 G k , then for a small enough, by Theorem 12.61 we have 
k =1 

/(a) = det (I n + aG(a)) 

(-l ) m+1 


= exp 


E 

m= 1 
oo 

ex p( E 


m= 1 
oo 


a m Tr(G(a) m )) 
m / 


-a m Tr 


ex P( E 


- 1 ) 


m 
m +1 


fc=l 


ra=l 


m 


oo 

E a kl+ '" +krn Tr(G kl ■ ■ ■ Gfe m )j) ■ 


(4.6) 


j”’ i^m — 1 


Since /(a) vanishes nowhere near 0, we can write /(a) = e 9 ^ a \ then by (14.61) . some direct compu¬ 
tation shows that 


<? (m) ( 0 )/m! = £ 
fe=l 


(_l)fc+i 

k 


E Tr(G Jt ---G h )). 

j 1-1 - \-jk= m 
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(4.7) 









For a small enough, let g(a) be the function satisfying 


det(S , 7 Q (T) - A I n ) = det(M - A I n ) ■ exp(g(a)), 


(4.8) 


then the coefficients g^ k \0)/k\ could be determined by (14.71) . And we have the following theorem, 
which is the main result in this subsection. 

oo 

Theorem 4.2. Under the above assumption, let g(a) be the function in . Let g(a) = c m° m 

m= 1 

be its Taylor expansion. Then 


Cm. — 


E 

fc=l 


(- 1 ) 


k +1 


E Tr (° 

j lH-1 -jk=m 


31 ' Gj k 


(4.9) 


where G k = (M — A I n ) MMj., and M = S^q{T), 

M k = f D(ti) [ D{t 2 )■■■ [ D(t k )dt k ■ ■ ■ dt 2 dt\, k £ N. 
Jo Jo 


10 Jo 

We only list the first 4 terms 


<o 


(4.10) 


gW(0)=Tr(G 1 ), 


g^(0)/2 = Tr(G 2 )-- 2 Tr(G 2 1 ), 


g {3) ( 0)/3! = Tr(G 3 ) - Tr(G 3 G 2 ) + ^Tr(G?), 


ff (4) (0)/4! = Tr(G 4 ) - \Tr{G\ ) - Tr(GiG 3 ) + Tr{G\G 2 ) - jTr(G 4 ). 


By the definition of G k , 


Tr(Gi) = Tr{M l M(M - AI 2 n) _1 ) = Tr 


D(s)ds ■ M(M - AI 2ri ) 


-l 


Tr(G 2 ) = Tr(M 2 M(M - A I 2 n) _1 ) = Tr 



D(a)dads ■ M[M - A/an)” 1 ) • 


Generally, 


Tr(G£*) = Tr([ l' T D( tl ) P D(t 2 ) ■ ■ ■ P" 

VL J0 JO Jo 


D{t k )dt k ■ ■ ■ dt 2 dti ■ M(M 



and Tr(Gj 1 ■ ■ ■ Gj k ) could be given similarly. 

In the case of Hamiltonian, there are some symmetry property of M k . please refer [3] for the 
detail. 
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4.2 Trace formula for linearly parameterized first order ODE 

For D a = Do + aD, let D a take place of D in both sides of Hill-type formula (11.41) . we get 


det 


d \ ( d 

--D a + ul n \ [ Jt +P* 


= (-l) n \C(S)\e-^e~2fo T P D P dt det{S lDa {T) - e vT I n ). 


(4.11) 


As we have proved that, both sides of (14.111) are analytic functions on a. Notice that the left hand 
side 


±-D a + vP (j t + P« 


-i 


— ( — Do + aD + vl n j f — — Do + vl n 


I-Do + Vln) (j t +P0 


-1 


(4.12) 


Hence 


det 


dt 


D a + vl n 


dt 


+ Pc 


-l 


= det 


dt 


— Dq — aD + vl, 


dt 


— Dq + vl n 


-l 


• det 


!L-D 0 + V D (4 + p 0 


-l 


= det I id — aD ( — — Do + vln 


d 


-i N 


• det 


4- Do + ,/„) (| + p 0 


-1 


(4.13) 


Let 


F — D ( — — Dq + vl n 


-l 


then the left hand side of (14.lip 
d 


dt 


D a T vl n 


f(a) = det 
= det 
= det (id — aF ) • det 


dt 


+ Po 


dt Da + Uln )(dt Do + uIn 


-l 


det 


“77 — Do + vl n 

dt 


Tt + P ° 


-1 


— Dq + vl n 


dt 


+ Po 


(4.14) 


Notice that F is a Hilbert-Schmidt operator with trace finite condition with respect to Pk- By 
Theorem 12.61 for a small, 


-1 


f{a) =exp V — a m Tr{F 

\ z ' m 


det 


\m=l 


— Do + vl n 


+ Pc 


-1 


(4.15) 
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On the other hand, by Theorem 14.21 the right hand side of (14.1111 equals to 

(—l) n \C(S)\e~ !1J ^~e~^ T r (D a )dt det(5 , 7 Q! (T) - e uT I n ) 

= (—l) n \C(S)\e~ :ii £- e _ 3 fo Tr ( D o) dt det (M — XI n ) e~% fo Tr ( D l dt e 9 ^ a \ 

OO 

where g(a) = c m a rn satisfies 

m =1 

det(5 7Q (T) - e vT I n ) = det (M - A I n )e 9 ^\ 
and c m are given in ()4.9I) . Comparing (14.151) with (14.161) . we have 

9(a) = f [ T Tr(D)dt+Y J —a m Tr(F m ), 

Z lr\ Vfl 

• /u m= 1 

and hence 

Tr(F) = \ £ Tr ( D (t))dt - Ci, 

and 


(4.16) 


(4.17) 


Tr(F m ) = —mc m , m > 2. 


(4.18) 


Thus we get 


Theorem 4.3. Let i' € C such that ^ — Do — v is invertible, F = D (-^ — Dq + ul n ) \ f/ien 


rr(F) = i j\r(D(t))dt - Tr(Gi), 
and for any positive integer m > 2, 


Tr(F m ) = 


(- 1 )' 


k =1 


£ Tr[G h ---G jk ) 

jH- \-jk= r m 


(4.19) 


(4.20) 


where G & = (M — A^n) and M = £70 (T), 


r*i 


— 1 


M k = / ZA(ii) / D(i 2 ) ■ ■ ■/ D(tk)dtk ■ ■ ■ dt 2 dt 1 ,k € N. 


(4.21) 


For large m, the right hand side of (14.20[) is a little complicated. However, for m = 1,2, we can 
write it down more precisely. 


Corollary 4.4. Under the assumption as in Theorem\4-.3, 


Tr(F) = J Tr(D(t))dt — Tr(^ J -f 0 1 (t)D(t)^ 0 (t)dt ■ M(M - e uT I n ) 


and 


Tr(F 2 ) = - 2 Tr^J 7 0 1 ( t ) D ( t ho(t) J 7 0 1 (s)D(s)'yo(s)dsdt ■ M(M 


-Mr, 


\-l 


+Tr 


o 

T 

7 0 " 1 (t)T>(t) 7 o(t)dt • M(M - e vT I n )~^ 


i o 
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Especially 


Tr\ D[- + vI r 


= Tr[ j r ^ D(t)dt ■ — S(S — e uT In) -1 


1 , 


~^ Tr (J o D(t)dt ■ (5 + e l ' T ){S — e uT In)- 1 )^ . (4.22) 


Taking derivative on both sides of (14.221) . we get 


Tr\D[j t +uI n 


—2 \ r T 

uT. 


= -Te ul Tr{ / D(t)dt ■ 5(5 - e uT I n )~ 2 


(4.23) 


5 Examples 

Krein considered second order system in [B] 

y" + A R(t)y = 0, y(0) + y(T) = y'( 0) + y\T) = 0, (5.1) 

where R{t) is continuous path of real symmetric matrices on M n . Set R ave = ^ / Q T R(t)dt and 


X(t) — / (R(s) — R ave )ds + C, 


(5.2) 


where C is a constant matrix which is chosen such that X ave = 0. Let A j be the eigenvalues of 
(15.11) . Krein get 


and 


= Tr{m)dt ’ 

E = | £Tr(X 2 (t))dt + ^Tr[(j\(t)dt) 2 ). 


(5.3) 


(5.4) 


In this section, we will give a generalization of Krein’s trace formula from our viewpoint. To 
simplify the notation, let A(v) = — (^ + v) 2 and denote by 


Uni!# 1 — 


R(t)dt , 


which is a constant matrix. From (14.231) . we have 

Tr(RA^)- 1 ) = -uT 2 ■ Tr{R ave • 5(5 - u)~ 2 ), 

By taking derivative with respect to v on both sides of (15. 5p . we get 

wT 4 

Tr(RA(y )~ 2 ) = Tr{R ave S(S 2 + 4w5 + w 2 )(5 - w)" 4 ). 

b 


Especially, if / 0 T Rdt = 0, then 


Tr(RA{v)~ 2 ) = 0 


(5.5) 


(5.6) 


(5.7) 
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Please note that (|5.5D is obtained by taking derivative on both sides of (14.221) . Obviously, in this 


case, we could calculate Tr^D^jjfc + vln'j ^ for any k G N by taking derivative with respect to 
for both sides of (14.221) . 

[rA{u)~ 1 

(rA{u)~ 1 


The purpose of the remaining part of this subsection is to compute Tr ( RA{v) 1 ) .In fact, 

2i 

as the form of multiple integral. 


we will use our trace formula (14.201) to express Tr 
Inspired by Krein [6j, let 

X(t)= [ ( R{s) - R ave )ds + C , 

Jo 

where C is a constant matrix. Obviously, 

X(0) = X(T) = C. 

At first, we will calculate Tr R — R ave )A{v)~ l ^j 

Proposition 5.1. For any constant matrix C which satisfies CS = SC, X(t) is defined as above, 
then 

2n 


Tr 


((R - Rave)A(v) 

L 


= —2 wTr 

—4 Tr 


X 2 dt • S(S - u) 

t 


-2 


+ 2 Tr 


X(t)dt -S(S-uj)- 1 ' 2 


o 


[ X(t) [ X(s)dsdt- S(S -co)- 1 

Jo Jo 


(5.8) 


Proof. Recall that A(o) = — + v) 2 . Please note that for any bounded operators A, B such that 

AB is trace class operator, we have that Tr(AB) = Tr(BA). Hence 


Tr 


R - R ave )A(v) 


-i 


= Tr 


= Tr 


(R ~ Rave) + v'j {R— Rave)(y^ + 


d \ “ 2 

U 




Noting that X(t) = f 0 (R(s) — R aV e)ds + C for some constant matrix C, we have 


dt 


+ v)X - X[ — + iA = X = R- R, 


thus 


XI 4 + X' - (i + X'x = (i + - RaAf. + d" 1 


dt 


dt 


dt 


dt 


(5.9) 


(5.10) 


(5.11) 


By (15.91) and (15. lip . 


Tr 


R Rave)A[y 


,-l 


2i 


= Tr 
= TTr 

= TTr 


Hi 


+ "> -IJJ + ..) x 
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1 

0 / d \ — 2 

- 2 Tr 



+ TTr [X 2 A{u) 


-n 


(5.12) 








































From (15.51) 


Tr 


X 2 A(v) 


-i 


= —tvTTr 


X 2 dt ■ S(S - u) 


-2 


(5.13) 


To continue, we should calculate Tr 


X 


wy 


Corollary 14.41 is 0, thus 7o(t) = Tn M = S. It follows t 
Tr 


by using Corollary 14.41 In this case, Dq in 
rat 


= Tr 




X{t)dt ■ S(S - u)- 1 ' 2 

o 


- 2 Tr 


[ X(t) [ X(s)dsdt ■ S(S - tu)- 1 

Jo Jo 


By substituting (|5.14p and (15.131) into (15.121) . we have the desired result. 
Please note that if S = ±/ n , then we have 


Tr 


[ X(t) [ X(s)dsdt ■ S(S 

Jo Jo 


- 0J 


\-l 


1 

= -Tr 
2 


X(t)dt) -S(S-uj)- 1 


. (5.14) 

□ 

(5.15) 


Moreover, in this case, the constant matrix C in Proposition 15.II could be chosen arbitrary. Partic¬ 
ularly, C could be chosen such that X(t)dt = 0. 


Corollary 5.2. In the case S = ±/ n? if we choose C such that j'^X(t)dt = 0, then we have 


Tr 


((R - Rave)A{v) 


-1 


2 oj 


(1^) : 


;Tr 


X 2 dt 


(5.16) 


Proposition 5.3. Suppose R aV eS = SR ave , then 

TrjRAin)- 1 ) 2 ) 

wT 4 


6 

+2Tr 


Tr(R z ave S(S 2 + 4 uS + cu 2 )(S - u )~ 4 ) -2 uT ■ Tr ( / X 2 dt ■ S(S - w) 


\—2 


\ ^ 

X(t)dt • S(S - u)- 1 1 


— 4Tr 


X(t) J X(s)ds ■ S{S - 17) 


Proof. 


( RA W 1 ) = ((R ave + (R-R ave ))A{v)- 1 ) 2 

= ({R - R aV e)A(v ) _1 ) + (^RaveA^)- 1 ^) 

+ (R - Rave)Rave+ RaveAfu)^ 1 {R ~ R aV e)A{uy 1 . (5.18) 

Please note that R aV eS = SR ave implies that Ayy 1 commutes with R aV e- We have 

TjRAyy 1 ) 2 ) = Trj(R- R ave )Ayyy 2 ) + Tr(R ave Ayy 1 ) 2 ) 

+2 Tr((R - Rave)Rave Ayy 2 ). (5.19) 
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Direct computation shows that 


Tr(R ave A(v)~ L ) 2 ) = Tr(R 2 ave A(v)~ 2 ) 

o;T 4 „ ,„ 2 c , c2 


6 


Tr(R 2 ave S(S 2 + 4 ooS + oj 2 )(S - w)" 4 ). 


Since Jq(R — R aV e)dt = 0, by (15.71) we have 


Tr 


(R Rave)RaveA(y ) 


-2 


= 0. 


(5.20) 

(5.21) 


Combining (j5.19D with (15.81) . (|5.20l) and (15.211) . the desired result is proved. □ 

Corollary 5.4. In the case S = ±I n , if we choose C such that / Q T X(t)dt = 0, then we have 


Tr((RA(v)~ 1 ) 2 ) 

±(1 ± 4cj + oj 2 )loT 2 

= —---- 1 r 

6(lTw) 4 



2 uT 

^ (1 Tw) 2 



(5.22) 


Remark 5.5. More specially, for the case that Krein considered, that is, let S = — I n and v = 0, 
in this case, uj = 1. By L5.5\) . we have Krein’s trace formula L5.S\) . Moreover, by L5.22\) . we have 


Ej| = f Trix ^ )dt+ l Tr l(l 

which is Krein s trace formula Qj. 
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